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Abstract 

We compute the non-equilibrium stress tensor induced by a heavy quark moving through weakly 
coupled QCD plasma at the speed of light and compare the result to = 4 Super Yang Mills theory 
at strong coupling. The QCD Boltzmann equation is reformulated as a Fokker-Planck equation in 
a leading log approximation which is used to compute the induced stress. The transition from non- 
equilibrium at short distances to equilibrium at large distances is analyzed with first and second 
order hydrodynamics. 

Even after accounting for the obvious differences in shear lengths, the strongly coupled theory 
is significantly better described by hydrodynamics at sub-asymptotic distances. We argue that 
this difference between the kinetic and AdS / CFT theories is related to the second order hydrody- 
namic coefficient Tjt. is numerically large in units of the shear length for theories based on the 
Boltzmann equation. 



I. INTRODUCTION 



The goal of the relativistic heavy ion programs at RHIC and the LHC is to create and 
study the properties of the Quark Gluon Plasma (QGP). Since the real time dynamics 
can not be directly probed with the Lattice QCD, the transport properties of QGP are of 
particular interest. There is a consensus that the experimental results on collective flow 
imply that the shear viscosity to entropy ratio of the Quark Gluon Plasma is remarkably 
small [T] 

Hr.'-^. (1.1) 
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Although it is difficult to reconcile this experimental ratio with a quasi-particle picture of 
the plasma, the success of resummed perturbation theory at describing lattice data on the 
pressure at somewhat higher temperatures suggests that a quasi-particle picture might pro- 
vide a qualitative guide to the plasma dynamics and form a basis for further approximation 
schemes 

However, in TV = 4 Super Yang Mills (SYM) theory at strong coupling (and large N^) 
gauge-gravity duality can be used to compute the shear to entropy ratio exactly [5l [6] , 

^ = -- (1.2) 

The tantalizing similarity between the experimental ratio and this celebrated theoretical 
result hints that a strong coupling limit (without quasi-particles) might provide a better 
starting point for understanding the plasma dynamics. At the very least, the strongly cou- 
pled M = A theory is an analytically tractable limit that provides a useful foil to perturbative 
calculations based on a quasi-particle description. 

The goal of this work is to compare the steady state response of non-abelian plasma at 
weak and strong coupling to an infinitely heavy quark probe moving at the speed of light. 
This is the simplest setup where the plasma response to an energetic probe can be analyzed 
in detail [7HT3]. At long distances, the non-equilibrium disturbance produced by the heavy 
quark probe thermalizes and forms a Mach cone and diffusion wake. The original motivation 
for investigating the Mach cone was the observation of an unusual structure in measured 
two particle correlations [El [T3]. Today, after the analysis of Alver and Roland (TB] and 
others [T71 - I20] . these unusual correlations are understood as the hydrodynamic response 
to fluctuations in the initial geometry and not as the medium response to an energetic 
probe. (The Mach cone picture also dramatically fails to explain current measurements in 
several ways - see for example Ref. [2lj and the conclusions of Ref. [22].) The goal of this 
manuscript is not to explain current measurements, but rather to examine the differences 
between weak and strong coupling, and to study the approach to hydrodynamics in both 
cases. Although the current paper has no immediate phenomenological goals, the medium 
response to energetic partons is currently being studied by all the experimental collaborations 
in various ways [2B]. Thus, this calculation, which analyzes the "jet" medium interaction 
precisely and determines a source for hydrodynamics through second order in the gradient 
expansion, may be useful for phenomenology in further studies. 

In the strongly coupled theory the stress tensor induced by a flnite velocity heavy quark 
was computed using the AdS / CFT correspondence O E] ■ The approach to hydrodynamics 
was analyzed as well as the short distance behavior [TOl [T2| |2^. In particular we will largely 



2 



follow (and to a certain extent extend) the hydrodynamic analysis of Ref. [12] to determine 
a hydrodynamic source through second order in the gradient expansion for the kinetic and 
strongly coupled theories. In the AdS/CFT calculation the lightlike f — )■ c limit was not 
analyzed due to various technical complications. (Here and below v is the velocity of the 



heavy quark.) As discussed in Section II B, it is possible to set f = c throughout the 



calculation by choosing a different set of gauge invariants. 

At weak coupling the hydrodynamic source has not been computed. Nevertheless the 
appropriate source for kinetic theory was determined in Ref. [13], and several estimates 
have been given for how this kinetic source is transformed through the relaxation process 
to hydrodynamics [23]. We have simplified the source for kinetic theory considerably and 
determined the plasma response at large distances by solving the linearized kinetic theory. 
After comparing the hydrodynamic solution at large distances to the full (leading-log) kinetic 
theory results, the appropriate source at each order in the hydrodynamic expansion can 
be computed. This part of the calculation employs a computer code developed by us to 
determine spectral functions at finite u and k [2S]- As a by-product of these spectral 
functions we determined the hydrodynamic transport coefficients that appear through second 
order in the gradient expansion in a leading log approximation. These parameters will be 
needed below to precisely determine the hydrodynamic source through second order. 

This work is limited to the analysis of the kinetics for a single heavy quark moving 
from past infinity. It would be quite interesting to follow the evolution of a parton shower 
initiated at time t = and the subsequent hydrodynamic response at late times. Although 
this transition has not been worked out, several of the most important ingredients have 
already been clarified [271429] . We hope to address the thermalization of full parton showers 
in future work. 



II. PRELIMINARIES 



We consider an infinitely heavy quark with f ~ 1 moving through a stationary high tem- 
perature plasma from past infinity. We will calculate the medium response in two model 
theories - pure glue QCD at asymptotically weak coupling and M = 4 SYM at asymptoti- 
cally strong coupling. Both theories are conformal in this limit and therefore the background 
stress tensor takes the characteristic form 

T^" = dia.g{e,V,V,V), with e = 3P. (2.1) 

The heavy quark moves in the z direction and imparts energy and momentum to the plasma, 
which ultimately induces a non-equilibrium response, 5T^^ . The non-equilibrium stress 
tensor (5T°° and 5T^^ are functions of cylindrical and comoving coordinates xt and xl 
where 

Xt = a/x^ + , = atan2(?/, x) , and xl = z — vt . (2.2) 

Rotational invariance around the z axis determines (T^^^T^v^ in terms of T^^^ 

T^%t,x) =T'^''^{xL, Xt) cos ^r, (2.3) 
T°^(t,a;) =T°^^(xL,XT)sin(^,. (2.4) 
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A. Kinetic Theory with a Heavy Quark Probe 



At weak coupling kinetic theory determines the response of the plasma to the heavy quark 
probe. To determine this response we linearize the Boltzmann equation for f{t,x,p) = 
Up + 5f{t,x,p) with Up = \l{epl'^° — 1) and will restrict the calculation to pure glue QCD 
in a leading log- approximation for simplicit}|^ The Boltzmann equation in this limit reads 



{dt ^vp-d^) Sfit, X, p) = C[f, p] + S{t, X, p) 



(2.5) 



where Vp = p and S(t,x,p) is the (to be discussed) source of non-equilibrium gluons 
produced by the heavy quark moving through the plasma. In a leading log(T/mD) ap- 
proximation, the linearized collision integral simplifies to a momentum diffusion equation 
supplemented by gain terms 



d 

C[f,p] = T^iA^ ( np{l + Up) 



dp 

g^CAml 



dp^ 
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log 



T 

rriD 



5f 



np{l + n. 



vJ J 



gain terms , 



(2.6) 



Here the Debye mass for a pure glue theory is 

d^p np{l + Up) 



m 



D 



2g'TA 



(27r) 



T 



(2.7) 



where Ta = Nc is the trace normalization of the adjoint representation. The diffusion 
equation should be solved with absorptive boundary conditions at p = so the number of 
particles is not conserved during the evolution |30]. Thus the microscopic theory encoded 
by this diffusion equation is conformal, and the only conserved quantities are energy and 
momentum. 

The gain terms are responsible for energy and momentum conservation. Specifically, the 
energy and momentum that is transferred (per time, per degree of freedom, per volume) to 
the non-equilibrium excess 5/ by equilibrium bath is 



dE 
~dt 
dP 

~dt 



= -Tha 
= -Tha 



np{l + Up) 



d^p 



g np{l + Up] 



d_ 
dp 



np{l + Up) 



(2.8) 
(2.9) 



as can easily be found by integrating both sides of Eq. ( 2.5[ ) without the source. This energy 
and momentum transfer by the bath drives additional particles away from equilibrium and 
ultimately fixes the structure of the gain terms: 



gain terms 



1 

6^ 



1^ 

p2 Qp 



p np{l + Up 



dE 1 



— np(l + Up) 



dP 

~dt 



(2.10) 



Including quarks would only lead to minor changes to our results as can be seen from Fig. 4 of Ref. [26] 
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where for subsequent use we have defined 



= I 77^^p(l + np) = —. (2.11) 



With the gain terms it is easy to verify that energy and momentum are conserved. Pre- 
viously we analyzed the linear response of this system of equations and determined the 
hydrodynamic plasma parameters in terms of fiA [26] 

^ -0.4613 — , (2.12) 



e + P fiA 



7]/ sT 



=6.32. (2.13) 



Since the microscopic dynamics is conformal, linearized, and only conserves energy and 
momentum (and not particle number), is the only second order hydrodynamic coefficient 
that appears to this order. If there are additional conserved quantities and the dynamics 
is not conformal then there are a multitude of coefficients that appear at second order - 
see for example [311 132]. Further it must be emphasized that these second order transport 
coefficients are insufficient to describe the decay of initial transients (non-hydrodynamic 
modes) |32]. 

The shear viscosity naturally agrees with prior results [331 El]- The fact that t^, is 
somewhat large ~ 6 compared to the viscous length is a generic result of kinetic theory 
[35} [36] . Finally, we note that records the transverse momentum broadening of a bath 
particle due to the soft scatterings, and is related to the soft part of jet-quenching q parameter 
in a leading T /m£, approximation [37], gsoft/2 = 2TfiA- Thus, the leading-log limit provides 
a concrete relation between r]/s and q. 

We now analyze how a heavy quark disturbs this system in the same leading log approx- 
imation scheme. The leading log energy loss of the heavy quark was computed long ago by 
Braaten and Thoma [381 139], with the result that the momentum transferred to the medium 
per time {i.e. minus the drag-force) is 



dt \dt' dt 

where Uriv) is the drag coefficient in a leading log approximation 

^^^^-^logm. (2.16) 



Svr \Tno J 

In the last line we have taken the f — )■ 1 limit of relevance to this work. As discussed more 
completely below, we have implicitly taken the coupling to zero before taking this limit so 
that radiative energy loss can be neglected. 

The drag force arises as equilibrium gluons from the bath scatter off the heavy quark 
probe and are driven out of equilibrium by the scattering process. This scattering produces 
a source of non-equilibrium gluons located at the position of the quark, 

S{t,x,p) = S{p)6'\x-vt). 
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Appendix [A| analyzes this scattering process {g + Q g + Q) and determines the appropriate 
momentum space source, S{p). In the hmit v = 1 the source has a simple form involving 
only two spherical harmonics 



S{p) = fiF 



np{l + rip 



2dAi 



B 



2 

p 



[1 + 2np 
T 



T 



p ■ V 



for V 



(2.17) 



With this source it is straightforward to integrate 2dA J (f^P^ over Eq. (|2.5|) and to verify 



that the stress tensor satisfies 



dt 



6-^{x - vt) . 



(2.18) 



Our strategy to determine the energy-momentum tensor is the following. We take the 
Fourier transform with respect to x of the kinetic equation in Eq. ( 2.5[ ) 



+ iVp ■ k)Sf{uj, k,p) = C[5f,p] + 2nS{p)5{uj - v ■ k) 



(2.19) 



and solve the Boltzmann equation in Fourier space. The technology to do this is based on 
simple matrix inversion as was documented in Ref. [26]. Then we calculate the stress energy 
tensor in Fourier space using kinetic theory. By Fourier transforming the stress tensor back 
to coordinate space, we can determine the energy and momentum density distributions. 
Additional details about this procedure are given in Appendix [Bj 



B. AdS/CFT with a heavy quark probe 

To describe the response of the Af = 4 plasma to a heavy quark probe we will follow the 
notations and conventions of Refs. pTlHO] which should be referred to for all details. Briefly, 
a heavy quark is described in AdSs with a trailing string. The energy and momentum gained 
by the medium as the heavy quark traverses the plasma is again parameterized with the 
drag coefficient /xj?(u) 

f-(f.^)-^<'')("^^)- 

This coefficient is found by determining the energy and momentum flowing down the string 
into the black hole [iTtiiH] 

..M^f^. (2.21) 

As in the weakly coupled case, the stress tensor in the strongly coupled Af = 4 theory satisfles 



Eq. (2.18) with the energy-momentum transfer rates given by the corresponding strongly 
coupled formulas. The deposited energy and momentum leads ultimately to a hydrodynamic 
response in the strongly coupled theory. The linearized hydrodynamic parameters through 
second order are [5l |35l Hi] 

7] 1 



(2.22) 



e + P AttT ' 

4-21og(2) ~ 2.61. (2.23) 



r?/(e + P) 
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According to AdS/CFT duality [15], strongly coupled SYM plasma is dual to the 5d 
AdS-Schwarzschild geometry, which has the metric 



'2 r 



-f{u)dt^ + dx^ + 



du^ 



(2.24) 



Here u is the radial coordinate of the AdS geometry with -u = corresponding to the 
boundary, L is the AdS curvature radius, f{u) = 1 — u'^/u\ with Uh = I/ttT, and T is the 
(Hawking) temperature of the plasma and dual geometry. 

The addition of an infinitely massive fundamental quark to the SYM plasma is dual to 
the addition of a string to the AdS Schwarzschild geometry, with the string ending at m = 
[16] . The presence of the string perturbs the 5d geometry according to Einstein's equations 
and the near boundary behavior of the metric perturbation encodes the changes in the SYM 
stress tensor due to the presence of the quark j37|. 

In the large limit the 5d gravitational constant ~ is small. Consequently 

the back-reaction of the string on the geometry can be treated perturbatively by solving 
the string equations in the background metric and subsequently computing the metric per- 
turbations sourced by the string. Solving the string equations of motion leads to the well 
known trailing string profile |^ H3] 



•^string (^; ^) ^ 



t + 



tan 



u 



Uh 



1 Uh-U 

-log — 

2 Uh + u 



(2.25) 



This string profile describes a quark moving at constant velocity v and has the 5d stress 
tensor 



Oi 



''05 



ut 



t 



ViVjF, 
U^Vi „ 



t 



where 



00 



u V + uh^f 



X 



string J; 



(2.26a) 
(2.26b) 

(2.27) 



27rL3^r^w- 
and A is the 't Hooft coupling. 

The five dimensional stress tensor of the string perturbs the background geometry. The 
information contained in the linear metric perturbation sourced by the trailing string can 
be conveniently packaged into fields which are invariant under infinitesimal diffeomorphisms 
[nmo]. Defining 



G 



MN 



G 



(0) 
MN 



MN 



(2.28) 



where G^^^ is the background metric ( |2.24[ ) and Hmn is the perturbation, and introducing 
a space-time Fourier transform 

HMN{t,X,u) = 



duj d^q 



HMN{(^,q,u)e 



2ti (27r)3 

we find two convenient diffeomorphism invariant fields [T2 



(2.29) 
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lUJ 



lUJ 



(2.30) 
(2.31) 



Here sums over repeated indices are implied with i,j running from 1 to 3 and a running 
from 1 to 2, ' denotes differentiation with respect to u and 



Ci = — g X (i; X g) , €2 = — v x q . 

q± q± 



(2.32) 



The field Zq transforms as a scalar under rotations and the field Zi transforms as a vector 
under rot at ions 

The equations of motion for Zq and Zi are straightforward but tedious to derive from 
the linearized Einstein equations. They read 



where 



Bo 
So 



Z'o + A-oZ'q + BoZo — KgS'o , 

_ 24 + Aq'^u^ + 6/ + q^u^f - 30p 
uf {vfiq^ + 6 — 6/) ' 
gV(14-5/-gV) + 18(4-/-3/^) 

'/' -I- 



u^f {q^u^ + 6 — 6/) 

{qH'^-?>q'q^)U 



8 4(gV+6-6/) , 



3MgV 



(2.33) 

(2.34) 
(2.35) 
(2.36) 



and 
where 



Siu 8u [g2 (gV+6) - / (12g2_9/")] gq^u 
+ ^ 3/2 (g2«2 - 6/ + 6) " ~T~ 



Z'l + AiZ[ + Bi Zi = kISi , 
_ uf - 3/ 



^55 — 8iqHi5 . 



3/2 - u {uq^ + 3/') / + u^co^ 



Bi 



Si 



f 



t'oi + iuJti5 e\t 



(2.37) 

(2.38) 
(2.39) 
(2.40) 



We note that since the string stress tensor (2.26) only depends on time through the 
combination x — vt, when Fourier transformed, the string stress tensor is proportional to 

^ A complete set of gauge invariants also includes a field which transforms as a traceless symmetric tensor 
under rotations J2 ■ This tensor mode determines the spatial components of the SYM stress and is not 
necessary for our purposes. 
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2n6{ijj — v-q). Consequently the f ields Zg are also proportional to 271 5 {uj — v-q). More over, 
becau se the string stress tensor (2.26) is proportional to — v"^ and Eqs. (2.33) and 

(2.37) are linear, we define Zg = Zsj\/\ — and solve for Z^ in the f — )■ 1 limit. 



Under the assumption that the boundary geometry is fiat, near the boundary the fields 
Zg have the asymptotic expansions 



^3 + ^(4)^4 ^ 



(2.41) 



The cubic expansion coefficients zf^ determine the SYM energy density (5r°° and the SYM 
energy fiux bT^^ via 



1 



1 

2^7r^ 



K3) 



g2 



(2.42) 

(2.43) 



For a given momentum q, to determine the SYM energy density and energy fiux we 
solve Eqs. (2.33) and (2.37) using pseudospectral methods. At the boundary at m = we 



impose the boundary condition that the fields have asymptotics of the form given in the 



series expansions (2.41), which is tantamount to demanding that the boundary geometry is 
fiat. At the horizon ai u = Uh we impose the boundary condition of infalling waves. This is 
tantamount to demanding that ~ (m — m/j)"*'^"''/^ near the horizon. With the Zg known, 
we can then extract the expansion coefficients Zg and construct the SYM energy density 



and energy fiux from Eq. (2.42) and Eq. (2.43) 



III. COMPARING ADS/CFT AND KINETIC THEORY 



Using the formalism outlined in the previous section we compute the energy density 
and Poynting vector induced by the heavy quark in both kinetic theory and the AdS / CFT 
correspondence. To compare the AdS/CFT and the kinetic theory results we have measured 
all length-scales in units of the shear length 



4 



(e + V)cl 



(3.1) 



where is the squared sound speed and in practice the speed of light is set to unity. Lo 
is proportional to the mean free path in kinetic theory and equal to l/vrT for the = 4 
theory. At large distances, where ideal hydrodynamics is applicable, the amplitude of the 
disturbance is proportional to the strength of the energy loss. Thus, we divide the response 
by the corresponding drag coefficient iif{v) for each theory, Eq. (2.15) and Eq. (2.21). With 



these rescalings the two theories produce the same (rescaled) stress tensor at asymptotically 
large distances, but differ in their approach to the ideal hydrodynamic limit as we will 
analyze in Section IV Fig. [T] and Fig. |2] compare the non-equilibrium stress in the two cases. 
A complete discussion is reserved for the summary in Section |V| 
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FIG. 1. The energy density (in scaled units) times R = + x\ that is induced by a heavy 
quark probe in (a) weakly coupled QCD and (b) strongly coupled A/" = 4 SYM. Here Lq is the 
shear length and the ^J-riv) is the drag coefficient for each case (see text). 



In order to compare the stress tensor quantitatively we plot the energy density in con- 
centric circles of radius R around the head of the quark. Specifically, we define 



dER 



2rf sin 9r (5T°°(R 



(3.2) 
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L o 



FIG. 2. The magnitude of the Poynting vector |T'^'| (in scaled units) times R= yx^ + that 
is induced by a heavy quark probe in (a) weakly coupled QCD and (b) strongly coupled = 4 
SYM. Here Lq is the shear length and the fJ-riv) is the drag coefficient for each case (see text). 

where R = xtXt + x^z and the polar angle is measured from the direction of the quark z: 
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Similarly, the angular distribution of the energy flux is given by 



^ = 27ri?2sin^^/?*5T°^(R), 

d^R 

= 2rfsin^R [cos9r5T^'(R) + sm9RST°''(R)] . (3.4) 

Numerical results for the angular distributions of the energy density and flux at several 
scaled distances UK = R/Lo are shown in Fig. |3} 

There is a dramatic change in the AdS/CFT curves between 5^ = 1 and = 5, indi- 
cating a transition from hydrodynamic behavior to quantum dynamics at distances of order 
~ l/vrT. Since this quantum dynamics lies beyond the semi-classical Boltzmann approxima- 
tion, no transition is seen in the kinetic theory curves. It would be interesting to calculate 
the stress tensor in this region perturbatively to better understand the differences between 
the two theories for i? ~ 1/ irT. 

Let us pause to discuss the limitations of both calculations. The point of the current 
work is to compare the approach to hydrodynamics at inflnitely weak and inflnitely strong 
coupling. In both cases the coupling is taken to zero or inflnity before the limit, f — )■ 1. 
As we now discuss, this limits the lengths scales that can be meaningfully studied in both 
theories. 

In the kinetic theory calculation the resulting stress tensor is valid for distances, R ^ 
1 / {g'^T log g~^) . For distances shorter than 1 / {g'^T log g~^) the collisionless non-abelian 
Vlasov equations should be used to describe the medium response at weak coupling [IHl US] . 
However, for distances longer than l/{g'^Tlogg~^) the effect of the plasma dynamics is in- 
corporated into the polarization tensor of the soft collisional integrals between the heavy 
quark and the particles that make up the bath. For example, the drag coefficient computed 
by Braaten and Thoma includes an HTL propagator in the t-channel exchange that includes 
the plasma physics of the polarization tensor [50j. We have limited the evaluation of this 
polarization tensor to a leading log approximation. 

Further, the weakly coupled calculation is limited to modest 7. We have implicitly taken 
the coupling constant to zero before taking f — )■ 1 so that the radiative energy loss of the 
heavy quark can be neglected. For a small but finite coupling constant radiative energy 
loss is suppressed when the Lorentz gamma factor of the heavy quark is not too large 

Similarly, the AdS/CFT calculation is limited to comparatively large distances xl,xt ^ 



l/y/^nT. In Section II B we introduced a new set of helicity variables so that the medium 
response on distances much greater than 1 / ^t^T can be determined in the limit when v = \. 
For distances much less than 1 / y/jT^T the structure of stress tensor has been analyzed in 
detail jlOl EH [52]. The medium response is characterized by a transverse length scale of 
1/y/^TiT, and a corresponding longitudinal scale of l/7^/^7rT, where 7 is the Lorentz factor 
of the heavy quark. In Fig. |3] we are investigating distances of order xl,xt ~ I/ttT and 
the physics associated with these very short scales is not visible. Although the importance 
of the l/^y^nT scale was understood in the context of momentum fiuctuations [531 - I55] . 
it is instructive to see these scales reappear in the asymptotic expansion for the induced 
stress tensor at short distances [TOl |23]. Rewriting Eq. (137) of Ref. |2^ in terms of 7, and 
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FIG. 3. The angular distribution of (a) the energy density [dEji/d9Yi]/ /j-f and (b) the energy flux 
[dSpi,/d^R]//iF at distances 9^ = 1, 5, 10, 20, and 40 for the kinetic theory and gauge gravity duahty. 
Here Lq is the shear length and the iJ-piv) is the drag coefficient for each case (see text). 



expanding for 7 large with xl = 'JXl of order xt, we have 

TV {2xI + xlxI) 



1 



Gvr^ (xr + x'i 



+ 



24 



™2 



[X 



+ X^) 



5/2 



+ . . . 



(3.5) 



The first term is the leading term at short distances and is independent of temperature. The 
second term (which captures the first finite temperature correction) is subleading in inverse 
powers of distance, but enhanced by a power of 7. Comparing the magnitude of these terms 
we see that the first term will dominate provided 



(7Xi 



~l~ Xrp ^ 



y77rT 



(3.6) 



This constraint limits the validity (and utility) of the short distance expansion to rather 
short distances. 

In summary we are examining two extreme limits - infinitely weak and infinitely strong 
coupling. The disadvantage of this approach is that some of the marked differences at short 
distances between the Vlasov response of weakly coupled QCD and the AdS/CFT response 
are not visible (see especially [56]). The advantage of this approach is that the onset of 
hydrodynamics can be clearly compared. We will analyze the hydrodynamic limit in the 
next section. 



IV. HYDRODYNAMIC ANALYSIS 



At large distances (see Fig. [T] and Fig. [2]), the medium response to the heavy quark probe 
clearly exhibits hydrodynamic flow. Following in part the discussion by Chesler and Yaffe 
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|12j . we will analyze this hydro dynamic response order by order in the gradient expansion 
for kinetic theory and for the AdS/CFT correspondence. The strongly coupled A/" = 4 
theory is conformal and the appropriate hydrodynamic theory is conformal hydrodynamics 
[351 [36] . Similarly, to leading order in the coupling constant QCD is also conformal and 
again conformal hydrodynamics is applicable in this limit. Beyond leading order, there are 
corrections to kinetic theory which break scale invariance and non-conformal hydrodynamics 
must be used to characterize the long wavelength response [2S1 1^ - 

For both kinetic theory and the AdS/CFT, the stress tensor of the full theory satisfies 



the conservation law Eq. (2.18) At large distances the stress tensor is described by hydrody- 
namics up to uniformly small corrections suppressed by inverse powers of the distance. In 
hydrodynamics, the spatial components of the stress tensor are specified by the constituent 
relation order by order in the gradient expansion. Specifically, the stress tensor can be 
written 

^Cdro = (e + V)u^u^ + Vg^^ + t^^^ , (4.1) 

where the dissipative part of the stress tensor 'k^'^ is expanded in gradients of and or 
(T and u^) to a specified order. For linearized conformal hydrodynamics (where = (1, u)) 
this expansion through second order in gradients reads [35] 

tt'^ = -2r] (d'u^) - 2r]T^ (d'd^ In T) (Static) , (4.2) 

and all temporal components are zero. Here ( ■ ■ ■ ) denotes the symmetric and traceless 
spatial component of the bracketed tensor [35J, i.e. for linearized hydrodynamics we have 

{d'd^ In T) = (^d'd^ - ^S'^d^^ In T . (4.3) 



We will refer to the conservation laws together with the constituent relation (Eq. (4.2)) as 
the static form of second order hydrodynamics. Using the lowest order equations of motion 
(ideal hydrodynamics) and conformal symmetry, the second order term {d^d'^ In T) can be 
replaced by the time derivative of vr^*^ [35J 

TT^^' = -2r]{d'u^) - T^dtn'^ (Dynamic) . (4.4) 

This rewrite of the constituent relation can be interpreted as a dynamical equation for n'^'^, 
and is similar to the second order form of Israel and Stewart [571 EH] • We will refer to this 
equation of motion for tc'^'^, together with the conservation laws as the dynamic form of 
second order hydrodynamics. 

At long distances, the form of the stress-energy tensor is described by T^yd^o ^P to terms 
suppressed by inverse powers of the distance. We will express the full stress tensor as a 
hydrodynamic term, which is irregular in the limit of w, fc — )■ 0, plus a correction which we 
will verify is a regular function for w, fc — )■ 

= Ti^y^,o[T'',T'']+r'r (4.5) 

We have temporarily emphasized here that Thydro ^ functional of the densities T''°,T°* as 
specified by the constituent relation and the equation of state. Then the equation of motion 
in Fourier space becomes 

- 6T'^ + ik^ 5Tg,,, = Si^^^S^, k) , (4.6) 
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where 



o3 
hydro 



dp-' 

[u, k) = —^2tt6{ijJ 



V ■ k) ~ ikW^ . 



(4.7) 



Examining ^hydro? we see that —ik^r^^ acts as an additional source term for hydrodynam- 
ics. What makes this decomposition useful is that r*-' (in contrast to Tj^y^ro) is a regular 
function at small k and u. For the steady state problem we are considering, r*-' can be 
written with three functions proportional to the symmetric tensors consisting of v and k 

t'^{u, = 27r/ip5(w -v-k) [ {vV - lv^6'^) 0i(u;, 

+ {tv'y + ik'v^ - qvik^S'^) 02 (t^, k^) + {k'k^ - \k^5'^) 03 (w, A;^) ] , (4.8) 

where 0i, 02 and 03 are regular for fc — t- 0. The source can be expressed similarly 

Shydro = 2t[5{uj — v-k) [0i,(u;, k'^)v + 0fc(w, k'^)ik\ . (4.9) 

Since t"^^ is localized, we can expand it for small u and k. Using an obvious notation for the 
Taylor series. 



+ 



dm, 



-ioj) + 



^(2,0) 



-lu) +01 



ikf +0{k^ 



(4.10) 



we see that the full source for hydro through second order can be expressed in terms of three 



expansion coefficients. 



^(0,0) 



and 



and 



'S'hydro = 2'Kllp6{uj — V-k) ^1 — 2(^01 



dm. 



iS^'V + 0r)A;2U 



=4>v 

K,2i(0,0) 1 ,2 ,{1,0) 



'I 



V 



'1 



lUJ 



iAm„ 

3</>2 



iuj I ik 



Oik^). (4.11) 



Summarizing, r^^ can be determined by comparing the full numeric al so lution for T*-' to 



Tu{j„„. Then, by fitting the functional form given by an expanded Eq. (|4.8l), we can extract 



■ hydro 

the three coefficients 



0i^'°\ and 02°''^'' for the Boltzmann equation and the AdS/CFT 
correspondence. These coefficients fully specify the hydrodynamic source of a heavy quark 
through quadratic order. Appendix [C] gives some sample fits to our numerical results and 
the fit coefficients are collected in Table IV The quality of the fits given in Appendix [C] 
indicates that r*-^ is well described by a polynomial at small k and u and justifies the analysis 
of this section. 



Examining Table |IV[ we notice that in the Boltzmann case the expansion coefficients 
proportional to 0i vanish. In fact, 0i(ci;, k"^) vanishes to all orders in u, k. This follows from 
rotational symmetry around the k axis and the somewhat special form of the kinetic theory 



source in Eq. (2.17). Since we do not expect this property to hold beyond the leading log 



approximation, further discussion of this point is relegated to Appendix [C} 

With the source functions 0t,(a;, /c^) and 0fc(a;, k"^) known (numerically) through quadratic 
order, the hydrodynamic approximation to the equations of motion for (static) second order 
hydrodynamics reads (with v = 1) 



-iu 5T°^' + c\k) ikST''^ + Tsk'' = [cos 6 (p^ + ik 0^] 27rfiFS{uj -v-k) 



-iu 5T°^'' + Dk'' 6T'>''' = sin 9 0^ 2n fipSiuj -v-k) 



nOx' 



(4.12) 
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Boltzmann 








0.484 


AdS/CFT 


-1 


-0.34 


-0.33 



TABLE I. Table of hydro dynamic source coefficients. The equations of motion are given by second 
order hydrodynamics with a so urce term, Eq. (4.6). The source term is expanded to quadratic 



order in k and uj in Eq. (4.11) which defines these coefficients. The first coefficient 



,(0,0) 



was 



computed analyticahy in the AdS / CFT case by Chesler and Yaffe pjj • Here Lq is the shear length 
(see text). 



where z' points along the k axis and x' is perpendicular to k (see Appendix |B]). In these 
equations = {Ar]/3)/{e + V), D = r]/{e + V). and c^{k) = cf(l + r^r^fc^). Using these 
approximate expressions and the exact equation 



27r j2p S{uj — V ■ k) 
the hydrodynamic solutions in Fourier space are given by 

+ k cos 9 (py] — k^ [Ts + 



5T°^''(a;,k) 
5T°^'(a;,k) 



i sin 9 1 



2n fip6{uj — k cos i 



271 hf5{uj — k cos I 



u + iD{u) P 
i [ UJ cos 9 (f)^ + c^{k) k] 



c^{k) k"^ -|- iVs ujk'^ 



2tt^fS{^ — k cos I 



(4.13) 

(4.14a) 
(4.14b) 
(4.14c) 



The solutions can also be used for first order hydrodynamics provided the wave speed c^(/c) 
and the source functions (p^ and (p^ are truncated at leading order, i.e. c^i^k) — )■ and 
01(0;, fc^) ~ (p^i'^K Similarly, the hydrodynamic solutions for the dynamic implementation of 



second order hydrodynamics takes the same functional form as Eq. (4.14) with the replace- 
ments 

r, „ . D 



c\k) ^ ci 



D Diuj) 



(4.15) 



1 — 'ir^cj 1 — iT-j^uj 

Given these hydrodynamic solutions and the hydrodynamic source functions tabulated 



in Table [TV] the hydrodynamic stress tensor in coordinate space can be computed using 
numerical Fourier transforms. The stress tensor of first and second order hydrodynamics 
(with the corresponding source) is compared to the full kinetic theory stress tensor in Fig. |4| 
Fig. [5] presents the analogous AdS / CFT results. Finally, a comparison between the static 
and dynamic implementations of second order hydrodynamics is given in Figure [6] and 
provides an estimate of higher order terms in the hydrodynamic expansion. We will discuss 
these results in the next section. 



V. SUMMARY AND DISCUSSION 



To keep this discussion self contained, we first recapitulate the problem and the corre- 
sponding notation. An infinitely heavy quark moves along the z axis with velocity f ~ c, 
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FIG. 4. The angular distribution of (a) the energy density [dER/d6'Fj,]/;UF and (b) the energy flux 
[dSpi,/d0R]//iF given by the Boltzmann equation at distances d\ = 10, 20, and 40. The Boltzmann 
results are compared with the first order and second order static hydrodynamics. 
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FIG. 5. The angular distribution of (a) the energy density did^\l [ipiv] and (b) the energy 

flux \diS n\l [ipiy) given by the AdS/CFT correspondence at distances 91 = 5, 10, and 20. The 
AdS/CFT results are compared with the first order and second order static hydrodynamics. Here 
Lo = l/vrT is the shear length and (J^riv) = 7\/A'?rT^/2 is the heavy quark drag coefficient for the 
AdS/CFT. 



depositing energy and disturbing the surrounding equilibrium plasma. We presented and 
compared the energy and momentum density distributions in two distinctly different plas- 
mas - a weakly coupled QCD plasma described by kinetic theory and a strongly coupled 
A/" = 4 plasma described by the AdS/CFT correspondence. The steady state stress tensor 



distributions can be written in cylindrical coordinates and comoving coordinates (see Eqs 2.2 
and 2.3). Fig. HI and Fig. exhibit the energy density and the magnitude of the Poynt- 
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FIG. 6. The angular distribution of the energy density [dER/d0R]//iF in (a) kinetic theory and (b) 
the AdS/CFT. The kinetic theory curves are plotted at distances 9^ = 20, 40 while the AdS/CFT 
curves are plotted at distances 9^ = 10, 20. The Boltzmann and AdS/CFT results are compared to 
the static and dynamic implementations of second order hydrodynamics (see text). The differences 
between the static and dynamic implementations of second order hydrodynamics reflects the size 
of neglected third order terms. Here Lq is the shear length and Hf{v) is the drag coefficient for 
each theory. 



ing vector |S| = |T°*| for kinetic theory and the A/" = 4 theory respectively. To compare 
these theories we measured all distances in terms of a length scale given by a combination 
of hydrodynamic parameters 

° (e + V)cj • 

This length is of order the mean free path in kinetic theory and equals l/vrT for the 
AdS/CFT. In each theory we divided the stress tensor by the corresponding heavy quark 
drag coefficient ^f{v) so that at asymptotically large distances (where ideal hydrodynamics 
is valid) the rescaled stress tensors are equal. At asymptotic distances both theories re- 
produce the "Mach cone" structure characteristic of ideal hydrodynamics, but these model 
plasmas differ at sub-asymptotic distances in their approach to this ideal hydrodynamic 
regime. In particular the Boltzmann theory is considerably less diffuse than the AdS/CFT. 
In kinetic theory the short distance response is reactive, and the sharp band seen in Fig. |2j^a) 
(which is indicative of free streaming quasi-particles) is absent in Fig. |2](b). We will see that 
the response of the AdS/CFT closely follows the predictions of hydrodynamics at modest 
distances which is strongly damped by the shear viscosity at short distances. This difference 
between the two can also be seen quantitatively in Fig. [3] which compares the kinetic theory 
and AdS/CFT results by plotting the energy density and energy flux at concentric circles 
of radius R in scaled units: 

= R/Lo . 
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The precise definitions of dEu/dOR and dSn/dOR is given by Eqs. 3.2, 3.3, and 3.4 



re- 



spectively. As discussed in the previous paragraph, the AdS/CFT curves are considerably 
broader than the corresponding kinetic theory results for 9i > 5. 

Examining Fig. |3] a striking feature of the AdS result is the dramatic transition from 
hydrodynamic behavior at = 5 to vacuum physics at = 1 which is not present in the 
kinetic theory calculation. This transition was noted previously and suggested as a way to 
reveal the strong coupling dynamics experimentally [52] . However, the absence of this transi- 
tion in the weak coupling calculation refiects a limitation of the kinetic theory approximation 
to QCD rather than a distinguishable difference between the AdS/CFT correspondence and 
weakly coupled QCD. Certainly the quantum dynamics at i? ~ I/ttT can not be captured 
by the semi-classical kinetic theory results. It would be interesting to compute the stress 
tensor in this region in fixed order finite temperature perturbation theory to see if the dy- 



namics of the two theories is similar at these length scales. As discussed in Section III using 



a short distance expansion [101 [21], in the AdS/CFT calculation new length scales emerge 
at distances of order l/^^/'^-wT and l/7^/^7rT which are not visible in Fig. [3} These scales 
have been associated with saturation physics [551 159] . 

Finally, we have analyzed the transition to the hydrodynamic regime in kinetic theory 
and the AdS/CFT. In particular we have determined the source appropriate for first and 
second order hydrodynamics in each theory, following in part a method outlined by Chesler 
and Yaffe [12]. This is elaborated in Section IV, and the precise source (which takes the form 



of derivatives of delta functions acting at the position of the quark) is given in Table IV 



The source is constructed so that the hydrodynamics to a given order, together with a 
source at the same order, reproduces the stress tensor of the full theory up to higher order 
powers of £mfp/i?. Fig. [s] given in Appendix [B] fits our extracted hydrodynamic source with 
a polynomial at small k and the superb agreement with our full numerical results justifies 



the source analysis of Section IV 



Fig. |4] compares first order and second order hydrodynamics to the kinetic theory results. 
Generally the second order theory provides only a minor improvement to the first order 
results until rather large radii, > 40. Indeed the behavior of the second order theory seems 
rather unphysical for 5^ < 10. This shows the limitations of second order hydrodynamics. 
Second order hydrodynamics is constructed to reproduce the full results order by order 
at asymptotic distances and is not constructed to describe the decay of non-equilibrium 
transients produced by the heavy quark. 

The slow convergence of hydrodynamic expansion to the full results of kinetic theory can 
be contrasted with the rapid convergence seen in the AdS/CFT results in Fig. |5} In the 
AdS / CFT case we see that first (second) order hydrodynamics describes the full result at the 
20% (4%) level for 9i ^ 5. The agreement with hydrodynamics is not as good as described 
earlier by Chesler and Yaffe. This is because we are describing a quark moving with velocity 
f = 1, and we have found that the deviation from equilibrium is noticeably larger than the 
V = 0.75 quarks studied by these authors. In addition the energy density distributions show 
larger deviations from first order hydrodynamics and were not studied previously. However, 
once (important) second order hydrodynamic corrections are included, the agreement with 
hydrodynamics is remarkable already at modest fH. 

We should mention that we have used the "static" version of second hydrodynamics which 
specifies tt'^'^ with a constituent relation analogous to the first order constituent relation [53] . 
Israel-Stewart type equations rewrite and interpret the constituent relation as a dynamic 
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equation using lower order equations of motion. This renders the system of equations hyper- 
bohc and causal, but mixes orders in the gradient expansion. We have compared the static 
and the dynamic theories for the kinetic and AdS theories in Fig. |6} Generally the Israel- 
Stewart type resummations do not lead to a significant improvement. Indeed, at smaller 
than shown in Fig. [6] Israel-Stewart type resummations can lead to spurious shocks which 
are not reproduced by the full result. The differences between the static and dynamic theo- 
ries gives an estimate of higher order terms, and this difference is smaller in AdS/CFT than 
in kinetic theory at the same 

Clearly, the convergence to the hydrodynamic limit is significantly faster in the J\f = 
4 theory relative to kinetic theory, even when lengths are measured in the scaled units 
described by 9^1. We remark that in the AdS/CFT the second order hydrodynamic parameter 
r,r is a factor of 2.5 smaller in scaled units than the corresponding kinetic theory parameter, 

=6.32, (Kinetic Theory) (5.1) 



r)/sT 

=4-2 log(2) ~ 2.61 . (AdS/CFT) (5.2) 

7]/ SI 

Based on these coefficients, it is natural to expect that the convergence to the hydrodynamic 
limit is faster for the Af = 4 theory than the corresponding kinetic theory. In theories based 
on quasi-particles and kinetic theory it is difficult to reduce the value of r^r in scaled units 
significantly [36]. Thus, it would seem that our principal result of this study is reasonably 
generic. Specifically, based on the model theories studied in this work we expect theories 
without quasiparticles to approach the hydrodynamic limit several times faster (in scaled 
units) than theories based on a quasiparticle description. From a practical perspective of 
applying hydrodynamics to various almost equilibrium phenomena of heavy ion physics [e.g. 
the hydrodynamic flow due to jets and other local disturbances) this factor of two can be 
quite important. 
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Appendix A: The kinetic theory source in a leading-log approximation 

The source of non-equilibrium gluons arises as gluons scatter off the heavy quark, g + Q ^ 
g + Q. The squared matrix element for this process is 



2dA 



16 



2{K ■ Pf M 



+ 



Q2 A{K-PY 



(Al) 



where K is the heavy quark momentum, P is the gluon momentum, Q = P' — P is the 
four momentum transferred to the gluon, and we have averaged over the colors and spins 
of the external gluon. In a leading log approximation only the (first) most singular term is 
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kept. The source of non-equilibrium gluons of momentum p is obtained from the Boltzmann 
colhsion integral for the g + Q g + Q process 

S{t,x,p) =S{p)6^{x - vt) 
r I A4|2 

= - I, J^^J^^(2vr)V(P,,t)[/p/.(l + Mil + - + + fk)] , 

(A2) 



where the heavy quark distribution = (27r)^5'^(fc — kH)5^{x — vt) is out of equilibrium. 
We expand the source in a spherical harmonic basis in the (x, y, z) coordinate system 



Sip) = ^Simip)Hi.mip; 



zx] 



Lm 



21 + 1 

An 



y^^Sio{p)Pi{cos9pk] 



(A3) 



and note that the Sim vanishes for non-zero m due to the azimuthal symmetry of the prob- 
lem. Using the orthogonality of (cos 6'pfc) and the phase-space parameterization and kine- 
matic approximations used to analyze the energy loss of heavy quarks [51j, the expansion 
coefficients can be written El 



Sioip) 



21 + 1 
47r 



dq 



du 

V 

-vq ^ 



(A4) 

where u is the energy transfer, q = p' — p the three momentum transfer, and is the 
azimuthal angle. The matrix elements in this parameterization are 



16p2(A;' 



0^2 



2dA 



2(1 — V cosOkp) 
(g2 - a;2)2 



(A5) 



where cos^^p is expressed in terms of the integration variables, u, q and [51]. Now we 
consistently expand out the integrand to quadratic order in u/T and q/T. This includes 
three types of terms: (1) an expansion of the distribution functions to quadratic order, 
(2) an expansion of the angle cos 6pk to linear order in q/T and (3) an expansion of the 
Legendre polynomial to linear order, Pi{x + 5x) Pi{x) + P'i{x) 6x. With the full expansion 
we explicitly integrate over the azimuthal angle (p aiid the energy u, observing empirically 
that all harmonics vanish for / > 1 when the velocity is lightlike. Further, the / = and 



/ = 1 harmonics can be done analytically leading to a simple answer recorded in Eq. (2.17) 

'g'CpNr 



Soo{p) 

Swip) 



1 



47rr2 
1 



Vl2^T^ [ 2dA 



2dA 

g'CpN, 



log 
log 



T \ I 2T 

— + 

rxiD) \ p 



1 + 2A 



mr, ) 



/p(l + /p)(l + 2/p) 



(A6a) 
(A6b) 



The leading-log simplifications described in this paragraph were observed previously when 
computing the shear viscosity. 



^ Specifically, we use Eq. (B21) of Ref. [5T]. However, we have interchanged the role of p and fc to be 
consistent with the notation used in this work. 
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Appendix B: Numerical Details about Kinetic theory and the Fourier Transform 



The goal of this appendix is to give some of the details of how the stress tensor is computed 
in kinetic theory. Most of the notation and strategy follows an appendix of Ref. p6] and 
this reference should be consulted for the full details. 

The linearized Boltzmann equation in Fourier space reads 

{-iu + ivp ■ k)6f{uj, fc, p) = C[6f, p] + 2'kS{p)6{uj - v ■ k) , (Bl) 

where Vp = p is the particle velocity and the vector k in the laboratory coordinate system 
is 

k = {k^,k'^,k^) = k{sm9kCos(pk, sin6'fc sin(yCfc, cos^fc)- (B2) 

In order to solve Eq. ( |B1[ ) numerically, it is convenient to introduce the Fourier coordinate 
system {x',y', z') where z' points along the Fourier momentum k 

x' =^k x{vxk), (B3) 
kr 

y'=^vxk, (B4) 

z =k . (B5) 

Following our previous work [20] we re-express the source and ultimately the solution 6f in 
terms of real spherical harmonics with respect to the x',y',z' coordinate system: 

{1 for m = 

y/2 cos rrupp for m > , (B6) 
\/2sm \m\(pp for m < 

where Nim is a normalization factor [26] . We note that the unit vector p has the following 
components 

p^' = ^J^Hnip; z'x') , py' = z'x') , p'' = ^io(p; z'x') . (B7) 

Since the distribution function 6f is independent of the azimuthal angle of k with respect 
to the original x, y, z coordinate system, we will choose this azimuthal angle (pk to be zero 
so that k lies in the x,z plane. Then in the {x',y',z') coordinate system the vector v has 
the components 

V = (f^',f^',t;^') = (sin 61, 0, cos6') , 

and 

vp= \l^(^cos9Hio{p]z'x') + sm9Hii{p-z'x')^ . (B8) 

The steady state solution to the linearized Boltzmann equation is also expanded in the 
spherical harmonics defined above 

6f{u,k,p) = ^2tt6{u - V ■ k)np{l + np)ximip, k)Hirnip; z'x') . (B9) 

Im 
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and the Boltzmann equation for xim becomes 



'-iuj5iv + ikCl^)p np{l + np)xvm = Cim[Sf,p] + V-p np{l + np)—Sim{,P 



(BIO) 



where the m index is not summed over. Here C]^ is a Clebsch Gordan coefficient 
Cim[^fiP\ is the colhsion integral in this basis, the normahzation coefficient is 



n 



T^dA/J-A 



(Bll) 



and the source is 



Sim{p, 0) 



-2T 



p 



+ 




1 + 2n) {5ii5mo COS + Sndmi sin i 



(B12) 



For each value of {k^, k^) (in units of ^ia/T) the linear equations are solved for Fim = 
Xim/'H. Due to rotational invariance of the collision operator around the k axis, the matrix 
equation does not mix harmonics with different magnetic quantum numbers, i. e. the collision 
operator is diagonal in m. Thus, the matrix equation in I, I' is solved for m = 1 and m = 
separately. Harmonics with m > 1 are not sourced by the motion of the quark in a leading 
log approximation. 

After solving for Sf{LJ, k,p), the energy and momentum excess due to the moving quark 
can be computed using the kinetic theory 



6T'^{uj,k) = 2dA I ^^p^6f{oo,k,p), 

where T^^{u, k) is proportional to 27r 6{uj — v ■ k) 

6T°^{u, k) = 27t6{uj - v ■ k)f^^{k„ kr) . 

The relationship between the {x,y,z) and the {x',y',z') coordinate system is 

6T'^%uj,k) = cos e6T^'''{uj,k) -sine 6T^''{uj,k) 
6T^y{uj,k) =0 

5T'''iuj,k) =sin0 5T°"'(w,k) + cos^(5T°"'(a;,k), 



(B13) 

(B14) 

(B15) 
(B16) 
(BIT) 



In presenting these formulas we have taken k in the x, z plane, i.e. = 0. More generally 
rotational invariance dictates that (T^^^T^y^ is proportional to k^ 



6f^%k) =5f^^^{K, kr) cos^k 
5f°^(fc) =5T°^'^(A;„fcT)sin(^fc, 



(B18) 
(B19) 



and the preceding discussion with ipk = suffices to determine T'^'''^{kz, ^t)- 

The stress tensor is tabulated in the k^, kx plane and then Fourier transforms can be used 
to compute the stress tensor in coordinate space 

d^fc Mk-x 5T°^(a;, k) . (B20) 



5T°^(t, 



X 



dtu 
2^ 



{2nY 
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Employing the familiar identity 

^ikrxTcosi^r-^,) ^ j^^krxr) + 2 J2 ^''UkTXT) cos(n(<^, - ^k)) , (B21) 

n 

it is not difficult to show that 

ST'%x,,xt)= -^UkTXT) ^e^'^^6f'\h,kr), (B22a) 
Jo 2tt 

rOO I, 17, poo II, 

6T'^nx„XT)= -^JiikTXT) "^e^'^Hdf'^-iKM). (B22b) 

ST'^{x,,xt)= -^UkTXT) ^e''^'6f'\h,kT). (B22c) 
Jo -^TT 27r 



The Fourier integrals in Eq. (B22) are not particularly easy. In order to get a convergent 
integral, we ffist multiply the numerical data by a window function which eliminates the 
contributions of high frequency modes. For kinetic theory results, a sample window function 
is 

erf ((A; - A;max)/cr) 



w{k) 2 



(B23) 

with kraax = 7-5^a/T and a = 3.5/iyi/T while for the AdS/CFT the fcmax and cr was consid- 



erably larger, fcmax = SOvrT and a = GOnT. We computed the integrals in Eq. (B22) using 
two methods. The ffist method used brute force summation to compute these integrals. In 
the second method we re-parameterized the stress tensor by the source functions 0^ and 0^ 
used in the hydrodynamic analysis. Specifically, we define the functions (pv{k) and (pk by 
the equations 

6T'^'{uj,k) ^ ^sing0 (fc) 2Tr^^p6{io-kcose), (B24) 
u + iDk'^ 

(c, k) = ,._,.(,)p^,r.c.fc^ 2./., 6iu - k cos e) , (B25) 

where c^(fc) = c^(l -|- TTrTsk"^) and T^^{uj,k) is determined by energy-momentum conserva- 
tion. This reparametrization is used for all values of the Fourier momentum k. Such a 
reparametrization is useful because, in contrast to T^^ and , the functions (j)v{k) and 
(f)k{k) are smooth and can be easily and accurately interpolated. Then the integrals are 
computed using Gaussian quadrature with break points at the hydrodynamic poles. The 
brute force summation and the more sophisticated numerical scheme give the same answer 
at the 1% level and are independent of the cutoff parameters. 



Appendix C: Details on the Hydrodynamic Source in Kinetic Theory and the 
AdS/CFT 

The purpose of this appendix is to explain in somewhat greater detail how the coefficients 



given in Table IV are computed for both kinetic theory and the AdS/CFT. In the process, we 



will exhibit several fits to our numerical results. The quality of these fits indicates that the 
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deviation of the stress tensor from its hydrodynamic form at small k and u is well described 
by a multivariate polynomial and justifies the hydrodynamic analysis of Section [TV| We will 
first discuss the AdS/CFT theory and then indicate how the analysis can be applied to 
kinetic theory. 



1. AdS/CFT 

The applicability of second order hydrodynamics to the AdS/CFT has been questioned 
based on the analytic structure of retarded stress tensor correlators in the lower half plane 
[60| l6T] . However, provided second order hydrodynamics is used to describe the behavior of 
hydrodynamic modes (z.e a pole in the retarded Green function arbitrarily close to the real 
axis), rather than to model the decay of non-hydrodynamic modes or transients {i.e. addi- 
tional analytic structure in the lower half plane), second order hydrodynamics is applicable 
to this strongly coupled theory. This appendix serves to clarify this points. 

After computing the exact stress tensor of the full AdS / CFT theory we define the func- 



tions 0„(fc) and 4>k{k) as in Eq. (B24). This would seem to be simply a reparametrization of 
the original numerical data on T^^' and T'^^' with two functions (f)v{k) and (f)k{k). However, 
as we will see, the functions (j)v{k) and (j)k{k) are analytic functions of k while the original 



fields have poles arbitrarily close to the real axis as fc — )■ 0, see Eq. (4.14). The purpose of 
first and second order hydrodynamics is to describe the location of these poles. First order 
hydrodynamics determines the pole location to first order in /c^mfp, but neglects terms of 
order (/c^mfp)^ which are captured by second order hydrodynamics^ It should be empha- 
sized that the pole shift is a consequence of modifying the ideal hydrodynamic equations 
of motion by powers of fcfmfp rather than modifying the source. Since, the ideal solution 
has a hydrodynamic pole aX u = Cgk modifying the equations of motion does not simply 
correct the solution by simple powers fc^mfp close to the pole. We will determine the source 
functions (j)v{k) and (j)k{k) using first and second order hydrodynamics. Specifically, we de- 



termine (j)v{k) and (f)k{k) using Eq. (B24) (with the same numerical data for the full stress 
tensor and ), but in the first order case we set the second order transport coefficient 
to zero in these equation^ Tj, ^ 0. 

Since the source functions (p^ and (pk are functions of k and u = kcosd, we can expand 
these functions in Fourier series 

(f)^{k, cosO) = (f)v-o{k) + 2(j)^.i{k) cos6 + 2(f).,.2{k) cos 26' + ... . (CI) 
7TT(j)k{k, cose) = (pkdk) + 2(j)k;i{k) cose + 2(f)k-2{k) cos2^ + . . . . (C2) 

In Fig. [7]we plot the terms of the Fourier series (j)v;nik) and (pkinik) and fit these functions 
with a simple power law, k°'. 

As we see from the fit, (p^ is well described by a quadratic polynomial at small k 

k \ ^ , f k \ , f k ^ ^ 



0„;o(A;) = l-O.1643( — 1 , 0,., = 0.5z( — 1 , = 0.0870 ( — ) . (C3) 



4 



^mfp should be taken as 1 /ttT in the strongly coupled theory. 



^ We note that (t)y is the same for first and second order hydrodynamics. Only (f>k is affected by non-zero 
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FIG. 7. Hydrodynamic fits to the AdS/CFT source functions = 1 + Acl)y and 'n'T(pk = — g + A^^ 
at small k (see the text surrounding Eq. ( |C1| )). The Fourier coefficients displayed in this figure 
have been multiplied by a factor indicated in parentheses to increase visibility, and are fit with 
the functional form C/c". The dashed lines and open symbols have non- integer fit-powers and lie 
beyond the description of hydrodynamics to the specified order, i.e. the fit is not expected to work, 
(a) The n = 0, n = 1, and n = 2 Fourier coefficients of (p^. (b) The n = 0, n = 1, and n = 2 
Fourier coefficients of (pk, when <j)k is extracted using first order hydrodynamics, (c) Same as (b), 
but for second order hydrodynamics. 



For first order hydrodynamics the terms of quadratic order can be neglected. Our numerical 
result for (p^-i is nicely consistent with the first order analytic result Ref. [TT].) 

Now we examine kcpk using first and second order hydrodynamic^ When first order hy- 
drodynamics is used the source function is not well described by a polynomial (see Fig. [7](b)). 
However, we see that k(j)k decreases faster than k (as k^'^"^ for n = 1), and therefore k(f)k 

^ We discuss k(j)k instead of simply (j>k since the source for hydrodynamics is (j)k{k)k. 
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can be neglected in a first order hydrodynamic analysis of the long wavelength response to 
the heavy quark. A concerned reader might worry that the resulting source function k(j)k 
is not well described by a quadratic polynomial, and incorrectly conclude a local source for 
hydrodynamics can not be constructed beyond linear order. However, when second order 
hydrodynamics is used to extract the source through quadratic order (as is appropriate), 
k(f)k is well described by the quadratic polynomial - see the hnear fit in Fig. ^c) for (f)k-i. 
Numerically we find 

1 k 
0fc;O = -3, = O.lli— , (C4) 

up to non- analytic terms that fall faster than k"^. (These non-analytic terms could be 
removed by pushing the hydrodynamic analysis to third order.) In summary, we see that 
using second order hydrodynamics neatly removes the non-analytic behavior seen in the first 
order source, and then the source for second order hydrodynamics is then well described by 
a quadratic polynomial. The coefficient which shifts the hydrodynamic pole is universal 
- it was determined through an analysis of two point functions [35j , and the same coefficient 
determines the long distance response to the disturbance produced by a heavy quark. By 
contrast the source functions 0^, and 0^ are not universal but depend on the particular way 
in which the heavy quark couples to hydrodynamic modes. Using the fits in Eq. (C3) and 
Eq. (C4) (and the relation between 0t,,0fc and 0i,02 given in Eq. (4.11)), we parametrize 



this source by three numbers to quadratic order which are given in Table IV 



2. Kinetic Theory 

As in the previous appendix, we define the functions 0„ and (pk from the exact stress 
tensor (Eq. ( B24[ )), and again expand these functions in Fourier series as in Eq. (CI), but. 



as is appropriate for the kinetic theory calculation, vrT^^ is replaced with (/i^/T)0fe . Then 
the Fourier coefficients are fit with a polynomial which are shown in Fig. [8j It was verified 
that the other Fourier coefficients of (p^ and k(f)k that are not shown decrease faster than 
k"^ and thus lie outside of the hydrodynamic analysis which is restricted to second order 
inclusive. Examining the fits shown in Fig. |8} we see that the parameterization 

/kT\ ^ 

0,;o =1 + 1.66(1) — , (C5) 
1 kT 

0fc;i = --1.66(l)i— , (C6) 

6 /iA 

describes our numerical data at small k and u = k cos 6 quite well. The fact that (p^fi and 
(pk;! have the same fit coefficient (up to a symmetry factor of 1/6) is a consequence of the 
hydrodynamic analysis of Section |IV[ Comparing the fit with the functional form giv en i n 
the text (Eq. ( |4.1l| )), we see that we have a non-zero (p^2'^^ (which is recorded in Table |iv[). 



but the coefficients (pf''^'' and (p^i'^'' seem to vanish. In fact (pi{uj,k'^) vanishes to all orders 
in /c as we will now show. 



To this end we can return to the analysis given in Section IV, For a given k, we expect 
that there is a non-zero component of T^^{u, k) which transforms as a spin two tensor under 
azimuthal rotations around the k axis. Examining the decomposition of r*-' = T*-' — T^y^ro 
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FIG. 8. (a) A polynomial fit to the source function cpy-o defined in Eq. (CI). The fit result is 
recorded in Eq. (C5). (b) Using the fit coefficient from (a), a prediction is made cpk-,! (see Eq. (C5)) 



using the analysis of Section IV 



into tensor structures (Eq. (4.8)), and noting that hydrodynamics does not yield such a 
spin two tensor, we see that this spin two component of T^^ determines ^v^v^ — 0i, 
since this is the only tensor structure of r'-' that has a spin two component under azimuthal 
rotations around k. Examining the source for the Boltzmann equation given in Eq. (B12) 



of Appendix [Bj we see that the specific form of the source does not excite the spin two (i.e. 
m = 2) components of the distribution function 6f. Thus, since the Boltzmann equation 
does not mix harmonics of different spin, the spin two component of the stress tensor vanishes 
and so does 0i. This approximate symmetry is specific to the simplified form of the source 
which arises in a leading-log approximation and is not expected to hold more generally. 
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